Effect of scattering and contacts on current and electrostatics in carbon nanotubes 
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We computationally study the electrostatic potential profile and current carrying capacity of 
carbon nanotubes as a function of length and diameter. Our study is based on solving the non 
equilibrium Green's function and Poisson equations self-consistently, including the effect of electron- 
phonon scattering. A transition from ballistic to diffusive regime of electron transport with increase 
of applied bias is manifested by qualitative changes in potential profiles, differential conductance and 
electric field in a nanotube. In the low bias ballistic limit, most of the applied voltage drop occurs 
near the contacts. In addition, the electric field at the tube center increases proportionally with 
diameter. In contrast, at high biases, most of the applied voltage drops across the nanotube, and 
the electric field at the tube center decreases with increase in diameter. We find that the differential 
conductance can increase or decrease with bias as a result of an interplay of nanotube length, 
diameter and a quality factor of the contacts. From an application view point, we find that the 
current carrying capacity of nanotubes increases with increase in diameter. Finally, we investigate 
the role of inner tubes in affecting the current carried by the outermost tube of a multiwalled 
nanotube. 
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G 

O 

o 



> 

On 
CO 

o 
in 
o 

-I— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

Metallic carbon nanotubes are near ideal conductors 
of current 0, 0, IE S @- While a single nanotube 
can be used as an interconnect in molecular electronics, 
nanotube arrays have shown promise as more conven- 
tional interconnects in conjunction with silicon technol- 
ogy 0, Q ■ A single nanotube has two modes that carry 
current at the Fermi energy, which yields a low bias con- 
ductance of Ae 2 /h and resistance of 6.5 fcfi. This corre- 
sponds to a current of 155 /zA at a bias of IV. Noting that 
a nanotube diameter can be as small as 5 A, it is easy 
to estimate that an array of metallic carbon nanotubes 
can carry current densities larger than 10 10 A/cm 2 . In 
fact, current densities approaching 10 9 A/cm 2 have been 
demonstrated 

The diameter of both metallic and semiconducting 
nanotubes can vary from 5 A to many tens of nanome- 
ters, with the electronic properties determined by the 
chiral angle. The bandgap of semiconducting nanotubes 
decreases inversely with diameter. As a result, a large 
diameter nanotube with a radius of 19 nm will have a 
bandgap of less than 2.5 kT at room temperature, mean- 
ing that large diameter semiconducting nanotubes carry 
non negligible current. Further, electrical contact can 
be made to many shells of large diameter multiwalled 
carbon nanotubes. Reference [j| demonstrated a resis- 
tance of nearly 500 Q in a multiwalled sample, which cor- 
responds to about twelve conducting shells. Therefore, 
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both small and large diameter nanotubes are promising 
as interconnects. Experiments on small diameter car- 
bon nanotubes show that the differential conductance 
decreases with applied bias, for voltages larger than 150 
mV 0, 0, E| ■ Reference showed that the conduc- 
tance decreases with increase in bias was caused by re- 
flection of electrons incident in crossing subbands due to 
scattering with zone boundary phonons. The non cross- 
ing subbands of small diameter nanotubes do not carry 
current due to their large band gap. In contrast, large 
diameter nanotubes experimentally show an increase in 
conductance with applied bias [J, Q, U4, LUI- Refer- 
ence Q suggested that as the diameter increases, elec- 
trons may tunnel into non crossing subbands, thus caus- 
ing an increase in differential conductance with applied 
bias. The main drawbacks of the calculation in refer- 
ence 01 was that the results depended on the assumed 
form of potential drop in metallic nanotubes and further, 
electron-phonon scattering was neglected. Recently, we 
performed self-consistent calculations |15| , which showed 
a dramatic increase in differential conductance of large 
diameter nanotubes with bias, in the ballistic limit. In 
the current work, we present results for current flow and 
potential profile in metallic nanotubes from a more com- 
prehensive model, which includes both charge self con- 
sistency and electron-phonon scattering. The potential 
and current- voltage characteristics are studied as a func- 
tion of nanotube diameter and length. We are primarily 
interested in short (~ 100 nm) rather than long nan- 
otubes, where the physics is more interesting and tech- 
nological applications are promising. The nature of the 
metallic contacts is also important. From an experimen- 
tal view point the contact between a metal and a nan- 
otube can either be an end-contact or side-contact. The 
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end-contact corresponds to only the nanotube tip elec- 
tronically interacting with the metal contact. In exper- 
iments, the end-contacts usually involves strong chemi- 
cal modification of the nanotube at the metal-nanotubc 
interface Also, reference Q found that end-contacts 
without sufficient chemical modification of the nanotube- 
metal interface have a large contact resistance. Due to 
the uncertainty of the contact bandstructure, modeling 
experimental end-contacts even remotely correct is diffi- 
cult. The side-contacts correspond to coupling between 
metal and nanotube atoms over many unit cells of the 
nanotube, and can be thought of as a nanotube buried 
inside a metal. Most experimental configurations cor- 
respond to side-contacts U 0] . An important feature 
of the side contact is that the coupling between atoms 
in the nanotube is much stronger than coupling between 
nanotube and metal atoms, which means that the band- 
structure of the contacts is similar to that of the nanotube 
between the contacts. In the side-contact geometry, elec- 
trons are predominantly injected from the metal into the 
nanotube buried in metal and then transmitted to the 
nanotube region between contacts. In fact, as proof of 
such a process, scaling of conductance with contact area 
has been observed in the side-contacted geometry by ref- 
erences 0, 0] . Modeling has also shown that the con- 
ductance in metallic zigzag nanotubes can be close to the 
theoretical maximum of 4e 2 / /i, when there is sufficient 
overlap with the contact [17|. In this work, we consider 
the metal-nanotube contact in both the limiting cases of 
side- and end-contacts. 

The outline of the paper is as follows. In section ITU 
we describe the formalism. The electrostatics at low bias 
is presented in section IIII Al electrostatics at high bias 
and current-voltage characteristics are presented in sec- 
tion HUH The role of inner shells in affecting the poten- 
tial profile of a current carrying outer shell is described 
in section IlII CI End-contacts that form both good and 
poor contacts are discussed in section All Dl We present 
our conclusions in section llVl 



II. 



FORMALISM 



In this paper we consider only zigzag carbon nan- 
otubes. The analysis for armchair nanotubes is similar. 
The general form of the Hamiltonian for electrons in a 
carbon nanotube can be written as: 



H 



/ j L x,x' c i,x c 3,x' 
i,j,x,x' 



(1) 



The sum is taken over all rings i,j in transport direction 
and all atoms located at x, x' in each ring. 

We make the following common approximations: i) 
only nearest neighbors are included; each atom in an sp 2 - 
coordinated carbon nanotube has three nearest neigh- 
bors, located a cc = 1.42 A away; ii) the bandstruc- 
ture consists of only 7r-orbital, with the hopping param- 



U = e p = 0. Such a tight-binding model is adequate 
to model transport properties in undeformed nanotubes. 
Within these approximations, only the following param- 
eters are non zero: 



Ul = U„,Vi 



X 

i,i — 1 
x,x' 

M+l 



(2) 



t Sx,x>, Vi 



2k 



where a = a cc \fi. In (N, 0) zigzag nanotubes, the wave 
vector in the circumferential direction is quantized as q = 
2Ttq/Na, q = 1, 2, ...N, creating eigenmodes in the energy 
spectrum. By doing a Fourier expansion of c\ x and Cj lX in 
g-space and using Eq. we obtain a decoupled electron 
Hamiltonian in the eigenmode space: 

H = Y,H q (3) 

i 

i 

where 



2t cos(y) = tt,Vi = 2k (5) 
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- l q 




A+l,i 
~ 0. 



= t a — i2, Vi = 2k 



The 1-D tight-binding Hamiltonian H q describes a chain 
with two sites per unit cell with on-site potential U and 
hopping parameters t\ and £2 (Fig^). For numerical so- 
lution, the spatial grid corresponds to the rings of the 
nanotube, separated by a/2 with a unit cell length of 
3a/2 = 2.13 A, which is half the unit cell length of a 
zigzag nanotube. The number of spatial gridpoints Ny 
is equal to the number of rings. 

The subband dispersion relations are given by 



E q (k) = ±\t \\l 



,3ak. ,07r, o,<Z7Tn i 1//2 ,„n 

4cos( — )co S (^)+4cos 2 (^)\ (6) 



Therefore, when N = 3fc, there are two subbands with 
zero bandgap the tube is metallic. In the rest of the pa- 
per we distinguish between metallic or crossing subbands 
(q = N/3 and 2N/3) and semiconducting or non crossing 
subbands. 

For each subband q we solve a system of transport 
equations [Tsj : 



AG R,q 

AG < ' > '« 



^ E )) GA ' q 



(7) 

(8) 



E* M (£). The self- 



eter t a 



-2.77 eV and the on-site potential 



where A = E - H q — V „ 

energies E^ <l> (i?) and E^j^ :> (i?) represent the effect 
of contacts and scattering. 

The contacts are assumed to be reflectionless reservoirs 
maintained at equilibrium, i.e. they have well defined 
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chemical potentials, equal to that of the metal leads: Vs 
in the source and Vd in the drain. Further, the nanotubc 
and metal are assumed to have the same workfunction. 
The contact self-energies T,f '<'> =T^'<'> due 
to the source contact are found using the surface 
Green's function of the leads, which is the solution of 
the following system of equations: 



(ai 
{a 2 



t 2 n S 

l 292,q 
l l9l,q 



J ph,l,q)9l,q — 1 



yR 
ph,2,q 



)92, q 



1 



(9) 



where the indices 1 and 2 stand for the two sites of the 
unit cell, a^a = E — U — Vs and T,^ h 1 2 are electron- 
phonon self-energies at the first two nodes near the 
source. The Green's function for the drain are solved for 
in a similar way, by making the substitutions if 2 — > t 21 , 
9\,2,q ~ * 9\,2,q an d taking T^ h 12 at the drain end. The 
expression for contact self-energies is then given by: 



^S,q 

yR 

D,q 



t 2 a s 

l 292,q 

t 2 a D 

H9l,q 



= -2i/s3m[E«_] 



J S,q 
D,q 



(10) 



S> ? = 2i(l - /s)9fm[Sg 
E> , = 2i(l - /u)9fm[Eg ig ] , 

where /s,d are the Fermi factors in the source and drain 
leads. 

The electron-phonon scattering is treated within self- 
consistent Born approximation. The electron-phonon 
self-energies are due to elastic (acoustic phonon) and in- 
elastic (optical and zone-boundary phonon) scattering: 



yR,<,> _ yR,<,> 
ph 7 q e/,q 



J in,q 



Ki,q( E ) = J^ ^' ^" 1 ' 



(11) 



(12) 



^ik,q(E) — ^ Din,v,q,q> 
v,q' 

[(n B {hu v ) + l)G<' q \E + hw v ) 
+n B (huj y )G < > q ' \E - hu u )] , (13) 

^in,q(E) = ^ Din,v,q,q' 
v,q' 

[(n B (^ l/ ) + l)G > '«'(S-^ l/ ) 
+n B {hLu„)G > > q ' (E + hu;„)} , (14) 
3m[E* = [£> M - Z< h . q ]/2i . (15) 



E' -E 



dE' (16) 



The matrix elements squared due to particular scatter- 
ing mechanisms are chosen so as to satisfy experimentally 
measured values of the mean free path Xf' m . Reference 



[ijj reported Xf w 1.6 /jm and A™ « 10 nm for a tube 
with a diameter of 1.8 nm, corresponding to a (24,0) nan- 
otube. Inelastic scattering is due to zone boundary and 
optical phonon modes with energies of 160 and 200 meV. 
Since, matrix elements squared scale inversely with the 
chirality index 20] , one obtains for an (N, 0) nanotube 



D, 



1 



off 



,24 



el.in 



2tt DOS(E f ) ■ r, 



el.in 



(*?)• ( 17 ) 



The decrease of a scattering rate in a single subband with 
the chirality index is saturated by the eventual increase 
of the intersubband scattering due to the larger number 
of subbands. Thus, with the increase of the diameter, 
the mean free path of a nanotube approaches that of 
graphite. 

Electron charge and current density rii and Ji at each 
node i are found from the following equations: 



eVs + WkT jjp 

<%?m=- (is) 



A p r -eVs + WkT 

tE 1 



dE 



eV D -10kT 



Gf,iUE)tl +hi ^-. (19) 



The lower limit of integration of Eq. I|18|) is determined 
by V m in = — 3t a — Vd — AS m , where AE rn is a renormal- 
ization energy related to the real part of electron-phonon 
self-energy. 

We model the electrostatics of the nanotube as a sys- 
tem of point charges between the two contacts located 
at y = ys = and y = y B = L. The "perfect con- 
tacts" are modeled as parallel semi-infinite three dimen- 
sional metal leads that are maintained at fixed source 
and drain potentials: Vs for y < ys and Vd for y > y B . 
So, while the self-energies due to contacts are identical 
to that of a semi-infinite nanotube, the role of electro- 
statics is included by image charges corresponding to a 
perfect metal. The electrostatic potential consists of a 
linear drop due to a uniform electric field created by the 
leads and the potential due to the charges on the tube 
and their images 

V = -eVs - e{V D - Vs)(y - ys)/(yD - ys) 

+ J2G(i,j)( nj -N) (20) 

3 

with the Green's function 



G(i,j) = 



4A??re 



E E 



s /(y i -y j + 2nL) 2 +p 2 K 



y J(y i +y ] +2nL) 2 +pl l 



(21) 



Here, p^; is the radial projection of the vector between 
atom k at ring i and atom I at ring j. The summation 
is performed over all atoms I at ring j for an arbitrary 
value of k. Maintaining the nanotube atoms buried in 
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the metal at a fixed potential is close to reality because 
of the screening properties of 3D metals. Within a few 
atomic layers from the metal surface, the potential should 
have approached the bulk values. While the variation in 
potential in these few atomic layers of the 3D metal is not 
captured in our model, our conclusions on the nanotube 
electrostatics should not be significantly affected. 

The calculations for one bias point involve two simul- 
taneous iterative processes: Born iterations for electron- 
phonon self-energies [Eqs. I|7ll6|l ] and Poisson iterations 
of [Eqs. (|7l20[l ] for the potential profile and charge distri- 
bution. Typically, 3 — 5 Born iterations were performed 
before updating a potential profile using Eq. I|20[) . The 
solution of Eqs. (|7I8|) employs a recursive algorithm [Isj . 
which scales linearly with the number of nodes. 

A. The importance of self-consistency and proper 
treatment of electron-phonon scattering 

A major computational burden of the approach dis- 
cussed above are the self-consistent procedure itself and 
also the Kramers-Kronig relation [Eq. i|16|) ]. In order 
to take the integral in Eq. (|16fl . at first, we have to 
solve Eqs. (1711 5[) over the whole band [— 3t a — Vh — 
AE rn ;+3t + AE rn ]. Due to the presence of inelastic 
scattering, the energy grid has to be uniform. Typically, 
Ne = 10 4 energy grid points are used in order to achieve 
a required precision in computing charge. At second, 
Eq. H16|l requires Ne 2 Ny multiplications for each sub- 
band, which makes it very time consuming. Such com- 
putational requirements pose a question on whether and 
what kind of sophistication is required in order to obtain 
an I-V characteristics. For a case of small diameter nan- 
otubes, we note that a contribution to current by cross- 
ing subbands in metallic nanotubes under a moderate 
bias does not depend significantly on the potential profile. 
The reason is the density of states of crossing subbands 
is nearly constant around the Fermi energy and therefore 
the transmission is insensitive to the changes in the po- 
tential profile. It is also clear that the renormalization 
of subbands due to scattering affects only band edges, 
but does not influence the density of states of crossing 
subbands near the Fermi energy. These two facts, allow 
us to conclude that neither self-consistency nor Kramers- 
Kronig relation [Eq. (|16[) ] is necessary to obtain an I-V 
characteristics in small diameter metallic nanotubes un- 
der a moderate bias. The criterium for this approxima- 
tion is that bias is lower than the renormalized bandgap 
of the first non crossing subband Enci, given by 

E NC i = 2(\t \\l - 2cos(~ - ^)\ - AE rn ) (22) 

- 2(\t \TrV3/N - AE rn ), for large N. 

Such an approximation, applicable e.g. to a (12,0) nan- 
otube under a bias smaller than IV, while giving incor- 
rect potential profiles would still result in a correct cur- 
rent with or without scattering. 



As will be seen later, in studying current through large 
diameter nanotubes it is important to take into account 
a current contribution by non crossing subbands. Tun- 
neling current by non crossing subbands has a lower 
threshold bias [Eq. (|22l) ] and depends exponentially on 
the slope of the potential near the edges of the tube. 
Therefore, a self-consistent solution for the shape of tun- 
neling barrier is must for large diameter nanotubes at all 
biases. This, in turn, necessitates the exact knowledge 
of electron charge and density of states. We now discuss 
how much the real part of electron-phonon retarded self- 
energy $te[Ep h q (E)] and thus the renormalization of the 
density of states affects the potential profile. In previous 
studies [2lJ SftefS^ q (E)\ was set to zero, which naturally 
alleviates computational requirements. In Fig[21we show 
the density of states of crossing subbands versus energy 
in a (12,0) nanotube under a zero bias. The Poisson it- 
erations were switched off and the potential V(y) = 0. 
Three curves represent different approximations to elec- 
tron transport: a ballistic case, when both real and imag- 
inary parts of electron-phonon self-energies are set to zero 
(dash-dotted line), a scattering case, when only imagi- 
nary part is taken into account (dashed line) , and a scat- 
tering case, when both imaginary and real part are non 
zero (solid line). For scattering cases, transport equa- 
tions were iterated for long enough to achieve conver- 
gence of electron-phonon self-energies, i.e. self-consistent 
Born approximation is satisfied. The area under the bal- 
listic curve Q = J DOS(E)dE is equal to 2 which is twice 
the charge per subband. The important consequence of 
taking into account only imaginary part of retarded self- 
energy is that the area Q decreases, resulting in ~ 5% loss 
of electron charge. When both real and imaginary parts 
are taken into account, electron charge is recovered due 
to the shift of the subband bottom. The loss of charge 
when the real part is neglected results in a completely 
incorrect potential profile when solving transport equa- 
tions self-consistently with Poisson equation [Eq. (|20|) ]. 
In FigEl we show potential profiles with and without a 
real part of electron-phonon self-energy when both Born 
and Poisson iterations have converged. The applied volt- 
age drops symmetrically across the nanotube, when the 
scattering is treated properly (solid line). When the real 
part is neglected, the profile shows a severe down shift 
due to a missing electron charge (dashed line). Such 
incorrect potential profile will result in large error in 
current due to non crossing subbands in large diameter 
nanotubes. Another source of error is an overestimated 
bandgap [Eq. (|22|) ] and higher threshold for the onset of 
current due to non crossing subbands. In this work, all 
our self-consistent calculations included Eq. (|16fl for the 
real part of electron-phonon self-energy. 
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III. RESULTS 
A. Electrostatics at low bias 

The mean free path due to scattering with acoustic 
phonons is in the range of a micron |9j, At low 

biases, i.e. biases lower than inelastic phonon energy, 
electron-phonon scattering does not play a significant role 
in determining the potential profile for wires of moderate 
length (less than few hundreds of nanometer). The po- 
tential profiles for (12,0) nanotubes of lengths 21.3 and 
213 nm are shown in Fig. 0] The edges of the nanotubc 
near the contact rapidly screen the applied bias / electric 
field. The potential drop is divided unequally between 
different parts of the nanotube, with 90% of the applied 
bias falling within 1 nm from the edges for both lengths. 

While the density of states (DOS) per unit length of 
metallic nanotubes is independent of diameter, the DOS 
per atom is inversely proportional to the diameter: 

DOS{E F )/atom = , (23) 

where N is the number of atoms in a ring of an (N,0) nan- 
otube. As a result, we find that the screening of metallic 
nanotubes degrades with diameter. The potential drop 
for two nanotubes with diameters of 0.94 nm [(12,0) nan- 
otube] and 18.8 nm [(240,0) nanotube] are shown in Fig. 
|5j Clearly, screening is poorer in the larger diameter nan- 
otube. In fact, while the potential drops by 45 mV in a 
distance of 1 nm from the edge for the (12,0) nanotube, 
the potential drop is only 17 mV for the (240,0) nan- 
otube. The inset of Fig. [S] shows a substantially larger 
electric field away from the edges of the large diameter 
nanotube. 

The electric field at the center of the nanotube as a 
function of length is shown in Fig. for the tube with a 
diameter of 0.94 nm. We find that for all diameters, the 
electric field decreases more rapidly than j,, where L is 
the length of the nanotube. The exact power law however 
depends on the diameter. If the computed electric field 
versus length is fit to -p-, the exponent a increases with 
increase in chirality. The value of a increases from 1.25 
to 1.75 as the diameter increases from 0.94 to 18.85 nm. 
Similarly, electric field versus diameter can be fit to D h , 
where b is in the range between and 1. 



B. Electrostatics at high bias and current-voltage 
characteristics 

At biases larger than 150 mV, the main scattering 
mechanism in defect free carbon nanotubes is electron- 
phonon interaction. Kane et aZ.jiJ found that emission 
of zone boundary and optical phonons arc the dominant 
scattering mechanisms. The potential profiles at high 
bias for 42.6 and 213 nm long (12,0) nanotubes are shown 
in Fig. [7| Due to the increased resistivity of the tube the 



potential profile drops almost uniformly across the en- 
tire nanotube, which is qualitatively different compared 
to the low bias (and no scattering) results of Fig. The 
potential drop at the edges of the nanotube accounts for 
only 30% of the applied bias. Because a bias drops in 
the bulk of the tube at the expense of the edges, the po- 
tential drop at the edges also decreases with increase of 
nanotube length, falling to 15% for the longer nanotube. 
Due to the diameter dependence of the scattering rates, 
the (12,0) tube is 20 times more resistive than (240,0), 
with mean free path of 5 nm for (12,0) versus 100 nm for 
(240,0). As a result, contrary to the low bias case, the 
potential profiles for larger diameter tubes show a larger 
voltage drop at the edges (inset of Fig. EJ ■ 

A major consequence of different potential profiles of 
large and small diameter nanotubes at low and high bias 
is the electric field at the nanotube center: the electric 
field is smaller in the small diameter nanotube in the bal- 
listic limit but the onset of electron-phonon scattering at 
high bias makes the electric field larger. This interesting 
reversal in electric field, demonstrated in Fig. [S]is ratio- 
nalized by noting that the potential drop in the (12,0) 
nanotube has become almost linear because the mean 
free path is much smaller than nanotube length, unlike 
in the (240,0) nanotube. 

So far electrostatics of nanotubes under applied bias 
was discussed. We now compare the potential profiles 
in the ballistic limit and with scattering, in a non zero 
electric field but at equilibrium. We choose an external 
electric field corresponding to a IV bias but the Fermi 
levels of the source and drain contacts are set to —0.5V. 
The electrostatic potential profiles in the two cases are 
almost identical as shown in FigEJ in sharp contrast to 
the non equilibrium case. The reason is, at equilibrium, 
following the Thomas-Fermi model, the potential profile 
should depend only on DOS at Fermi energy, which is 
unaffected by electron-phonon scattering. 

We now discuss the current voltage characteristics of 
small and large diameter nanotubes. For the (12,0) nan- 
otube, only two crossing subbands contribute to trans- 
port. The ballistic current increases linearly with applied 
bias and the differential conductance is Ae 2 /h. Scatter- 
ing by inelastic phonons causes current saturation and 
the decrease of the differential conductance (Fig. ITU)l . As 
the length of the nanotube is increased, so is the number 
of scattering events. The family of current - voltage char- 
acteristics show the transition from ballistic to diffusive 
transport regime. The current saturation at the value of 
25 [iA for the longest tube agrees well with experimental 
data. 

Current and differential conductance versus bias for a 
42.6 nm nanotubes and a wide range of diameters are 
shown in Fig^] At low bias, phonon scattering is weak, 
so the current and conductance are still close to the bal- 
listic limit. The saturation of current corresponds to the 
onset of inelastic phonon scattering and to the decrease 
of conductance at high bias. Differential conductance 
for small and moderate diameter nanotubes (A = 12, 
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36 and 60) is qualitatively similar because of the same 
number of conducting modes. Quantitatively, the con- 
ductance increases with diameter due to decreasing scat- 
tering rates. This transition from low to high bias regime 
is also present in large diameter nanotubes (N = 90, 120 
and 240). 

We note that current increases with diameter. In the 
ballistic limit, the self-consistently calculated current of 
(240,0) nanotube at IV is 310 /LtA and the differential 
conductance is almost 13e 2 /h. When scattering is in- 
cluded the current is decreased to 218.3 fiA which is 
much larger than the current carried by a (12,0) nan- 
otube of the same length, which is 45.7 fiA. In addition, 
the differential conductance versus bias exhibits a qual- 
itative change in shape. It is bell-shaped for small di- 
ameter nanotubes and transforms to U-shaped for large 
diameter nanotubes. As explained in [l5j . the increase 
of differential conductance with bias occurs due to in- 
jection of electrons from contacts into the low energy 
non crossing subbands. A schematics of Zener tunnel- 
ing process is shown in Fig ll2l when the bias becomes 
larger than twice the bandgap of the lowest non crossing 
subband, electrons can tunnel from valence band states 
(E < —Enci/2 + V(y)) in the source to the conduction 
band states (E > E NC i/2 + V(y)) in the middle of the 
tube (the channel) and also from valence band states in 
the channel to the conduction band states in the drain. 
The lowest non crossing subbands (q = A/3 — 1) in (90,0), 
(120,0) and (240,0) in the ballistic limit have bandgaps 
Enci = 331, 249 and 125 meV respectively and start to 
contribute to current at biases of twice these values. We 
note that the above picture is based on the assumption 
that the applied voltage drops symmetrically across the 
nanotube. In the case when one contact is much more 
resistive than the other and the applied voltage drops 
completely at one of the edges, the threshold for Zener 
tunneling is reduced, starting at biases equal to E^ci 
rather than 2Emci- Additionally, if the contacts are 
metallic (and not perfect nanotube contacts), then the 
threshold bias for tunneling into non crossing subbands 
will be further reduced by a factor of 2, an issue that will 
be discussed in section fill Dl In addition to Zener tun- 
neling between states with the same quantum number q, 
scattering also induces a phonon assisted tunneling from 
bonding states of non crossing subbands to crossing sub- 
bands. Although these processes require intersubband 
scattering, the threshold bias is twice as low. 

Another important feature of differential conductance 
occurs at zero bias in (240,0) nanotube (Tig lllfl . The 
zero bias conductance of (240,0) nanotube is larger than 
4e 2 //i because the non crossing subbands are partially 
filled and contribute to current: the first non crossing 
subband opens at an energy of 2.4 kT from the band 
center. This contribution is a simple intraband trans- 
port, determined by the population of the conduction 
band in the source and the valence band in the drain, 
but rather insensitive to the details of the potential pro- 
file. At slightly higher biases the non crossing subbands 



contribution to current saturates to a constant value and 
the contribution to differential conductance decreases to 
zero, while the total conductance decreases to 4e 2 //i. 

Increasing the length of a nanotube eventually makes 
the nanotube length larger than the mean free path. The 
current carried by a 213 nm long (240,0) nanotube is now 
86.6 [iA at IV. Note that the differential conductance of a 
long nanotube, shown in the inset of Fig ll3l has changed 
to bell-shaped, in contrast to the short nanotube case 
presented in Fig^] This is because when the nanotube 
length is many times the mean free path, the potential 
drop in the nanotube bulk becomes more linear and as a 
result the barrier width (Tig ll2|l for tunneling into non 
crossing subbands increases. The, non crossing subbands 
then do not carry significant current. Over all, the volt- 
age drop at the edges is a crucial factor determining the 
current contribution due to non crossing subbands. It 
is directly related to the height and width of the bar- 
rier for tunneling into non crossing subbands (FigEJ}, 
which depends on the scattering rate, nanotube length 
and screening capability. 

In the rest of the paper we discuss the role of electro- 
static coupling to inner shells of a large diameter nan- 
otube and the role of contacts. Specifically, we will see 
that both issues lead to increased tunneling into non 
crossing subbands because they lead to a decrease in bar- 
rier thickness. 



C. Electrostatic gating due to inner shells 

In multiwalled nanotubes (MWNT) current flows only 
through the outermost shell as the inner shells are not 
electronically coupled to the metal leads 0], and inter 
shell hopping is negligible |22j . Here, we are interested in 
the role of electrostatic coupling between the inner and 
outer shells, in affecting the potential profile and current 
carried by the outer shell. We find that the inner shells 
effectively decrease the barrier width for tunneling into 
non crossing subbands by providing additional screen- 
ing. Fig 1 131 shows the potential profile in the outermost 
(240,0) shell. The inner shells are chosen to be metallic 
nanotubes with chirality indices N = 225, 216, 210 and 
201, such that a separation between the walls is roughly 
3 A. Even when just one metallic shell (225,0) is added, 
the barrier width for tunneling into the first non crossing 
subband decreases dramatically from 2.97 to 0.97 nm. 
Adding more shells further reduces the barrier width by 
smaller amounts, and the effect saturates. The differ- 
ential conductance with the inner shells is qualitatively 
different as shown for the 213 nm long (240,0) nanotube 
in the inset of Fig^| We find that the differential con- 
ductance at higher biases is larger but does not show an 
increase with bias, when inner shells are included. 
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D. Influence of contact quality 

We have so far assumed that the contacts are per- 
fect, i.e. made of semi-infinite carbon nanotube leads. 
As mentioned in the introduction, this physically corre- 
sponds to carbon nanotubes weakly coupled to the metal 
in which they are buried. So, injection of electrons into 
the carbon nanotube lying between the contacts occurs 
from the carbon nanotube buried in the metal. As a re- 
sult of this, injection into non crossing subbands from 
the contact cannot occur at E=0. Note that selection 
rules don't permit injection of electrons from the cross- 
ing subband of the nanotube contact into non crossing 
subbands between the contacts. In this section, we re- 
lax this condition and consider injection from the source 
contact into all nanotube subbands. The drain contact 
continues to be a nanotube contact. To accomplish this 
in a phenomenological manner, we model the contact self 
energies in Eq. (fltl)) as, 

Eg, = -iapt 2 q , (24) 

where t q is the original hopping parameter in Eq. 1)100. 
a is called the quality factor of the contacts, and p is 
the density of states of the contact which is chosen to 
be close to that of gold, p = Q.lleV^ 1 . The self energy 
in Eq. (|24J) permits injection into the first non crossing 
subband at a bias of -E/vci rather than the 2£Vci with 
nanotube contacts. 

We now present results for the potential profiles with 
good (a = 1) and poor (a — 10~ 3 and a — 10 3 ) con- 
tacts in Fig^] Note that the terminology of poor and 
good contact is a relative to the intrinsic resistance of the 
nanotube that arises due to electron-phonon scattering. 
For a = 1, the potential profile of a (240,0) nanotube 
of length 213 nm is more or less symmetric and qualita- 
tively similar to the perfect contact results presented in 
Fig'ED For a — 10~ 3 and a — I0 3 , the electrostatic po- 
tential drops predominantly at the source-end due to the 
large source-end contact resistance. The large potential 
drop and hence extremely thin barrier in the source-end 
will facilitate tunneling into the non crossing subbands 
for poor contact. As mentioned before, the asymmetry 
of the potential profile further reduces the threshold for 
Zener tunneling to biases of Emci/^- The difference be- 
tween the potential profiles for good and poor contacts 
has a profound influence on the shape of the differential 
conductance versus bias as shown in Fig lI5l For good 
contacts, the 213 nm long nanotube shows a decreas- 
ing conductance with bias, in agreement with the per- 
fect contact case. In contrast to this, for poor contacts, 
the differential conductance increases with bias. This in- 
crease with bias is a direct consequence of the potential 
profile shown by the dashed line in Fig lI4l which facili- 
tates significant injection into non crossing subbands. To 
see more directly that the non crossing subbands are im- 
portant in carrying current in long nanotubes with poor 
contacts, we calculate the current carried at a drain bias 
of 0.5 V as a function of the quality factor of the contacts, 



with only crossing subbands and with all subbands. For 
good contacts, the current with only crossing subbands 
is very close to the current with all subbands (Fia lf (ill . 
But for poor contacts, the current with all subbands is 
significantly higher than the current with only crossing 
subbands. For example, when a = 10~ 3 , the current with 
all subbands is nearly an order of magnitude larger than 
with only crossing subbands. From Fig lltjl one can also 
see that the U-shaped differential conductance curve ver- 
sus bias should be observed in a wide range of a except 
a small region near a ~ I , where non crossing subbands 
do not carry a significant current. 



IV. CONCLUSION 

We have studied transport in nanotubes of varying 
diameters and length, with electron-phonon scattering 
included. We find that charge self-consistency and the 
proper treatment of subband renormalization due to scat- 
tering are crucial in determining the correct current- 
voltage characteristics in large diameter metallic nan- 
otubes. In the small bias ballistic limit, while the applied 
bias drops predominantly at the nanotube-contact inter- 
faces, screening is incomplete at the tube center. Further, 
screening improves with decrease in nanotube diameter 
due to the increased density of states per atom near the 
Fermi energy. At biases larger than 150 mV, electron- 
phonon scattering becomes important and the electro- 
static potential drops primarily in the bulk of the nan- 
otube rather than at nanotube-contact interfaces. How- 
ever, as the mean free path for electron-phonon scattering 
increases with increase in nanotube diameter, the poten- 
tial drop in the bulk of the nanotube is larger for small 
diameter nanotubes. As a result, the electric field at 
the nanotube center increases with increase in diameter 
at small biases, and decreases with increase in diameter 
at high biases. This interesting reversal in electric field 
versus diameter is computationally seen for nanotube 
lengths from 42 to 213 nm. Over all, we find that larger 
diameter nanotubes are capable of carrying more cur- 
rent because of increase in mean free path with increase 
in diameter and lower bias threshold for injection into 
non crossing subbands. For small diameter nanotubes of 
length 42 nm, we find that the differential conductance 
versus bias is bell-shaped, with the largest value at zero 
bias. The reason for the smaller differential conductance 
at larger bias is reflection of electrons due to inelastic 
phonon scattering and availability of only crossing sub- 
bands for charge transport. In contrast, large diameter 
nanotubes of the same length show an increase in differ- 
ential conductance with increase in bias due to injection 
into non crossing subbands and larger mean free paths. 
At a nanotube length (213 nm) longer than the mean free 
path, we find that the differential conductance for large 
diameter nanotubes transitions to a bell-shaped curve 
similar to the small diameter case, when coupling to con- 
tacts is good (small contact resistance). We have also 
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modeled the role of inner shells in affecting the poten- 
tial profile of the outer shell of a multiwalled nanotube. 
Here, we find that the inner shells cause a change in the 
electrostatic potential profile of the outer shell so as to 
make the potential drop sharper at the nanotube-contact 
interfaces. Our computational study has shown that the 
potential drop and differential conductance in metallic 
nanotubes is determined by an interesting interplay of 
diameter, mean free path, nature of contacts and bias 
threshold for tunneling into non crossing subbands. Fi- 
nally, the most important result of this paper is that a 
quality of contacts is a primary factor in determining 
the shape of the differential conductance versus bias in 
large diameter nanotubes. When the resistance at the 
nanotube-contact interface is larger than the intrinsic re- 
sistance of the nanotube, there is a large potential drop at 



the interface. This potential drop facilitates considerable 
tunneling into non crossing subbands and as result the 
differential conductance of the 213 nm long (240,0) nan- 
otube with poor contact increases with bias, in contrast 
to the case with good contacts. This finding gives a possi- 
ble explanation to the increase in differential conductance 
with bias seen in the recent experiment of reference |13| . 
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FIG. 1: Zigzag carbon nanotube and the corresponding ID 
chain. The hopping parameter between nearest neighbors in 
the nanotube is t . The 1-D chain has two sites per unit 
cell with on-site potential U and hopping parameters t\ = 
2i D cos(^) and t 2 = t a 

„x 10~ 3 

- - - 3(11=0,916=0,0=2.00 ;j 

3m*0,9te=0,Q=1.89 jj 

3m*0,9le*0,Q=2.00 \\ 




Energy (eV) 

FIG. 2: Density of states in a (12,0) nanotube under a zero 
bias, for a fixed potential V(y) = 0. The area under each 
curve Q = J DOS(E)dE gives twice the charge per subband. 
Including only the imaginary part of scattering self-energy 
(dashed) results in a loss of charge. 
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FIG. 3: Converged potential profiles versus position for (12,0) 
nanotube for the cases when the real part of electron-phonon 
self-energy is neglected (dashed) or taken into account (solid) . 
Neglecting the real part results in a severe down shift of po- 
tential profile. 
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FIG. 4: Potential versus position for (12,0) nanotubes of 
lengths 21.3 and 213 nm. The diameter of the (12,0) nan- 
otube is 0.94 nm and the applied bias is 100 mV. 



11 



Or 



-0.02 



(12,0) Diameter= 0.94 nm 

(240,0) Diameter^ 18.8 nm 

Bias = 0.1V 




50 



100 150 
Y (nm) 



200 



FIG. 5: Potential versus position for (12,0) and (240,0) nan- 
otubes, which have diameters of 0.94 and 18.8 nm respec- 
tively. The screening for the large diameter nanotube is sig- 
nificantly poorer. The inset magnifies the potential close to 
the nanotube-contact interface. 
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FIG. 6: This plot shows the electric field at the mid point 
versus nanotube length. A (12,0) nanotube at an applied 
bias of 100 mV is considered. The electric field decreases 
approximately as L~ 125 . 
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FIG. 7: The potential as a function of position is shown for 
(12,0) nanotubes of lengths 42.6 and 213 nm (solid). Shown 
for comparison are the potential profile in the ballistic case 
(dashed). Inset: Potential versus position for two different 
diameters. Note that contrary to the low bias case, the electric 
fields away from the edges is larger for low diameter tube. The 
length of the nanotubes are both 213 nm. 
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FIG. 8: Electric field as a function of diameter for nanotubes 
of lengths 42.6, 85.2, 127.8 and 213 nm. (a) At low bias, 
higher density of states in lower diameter tubes leads to the 
better screening (lower field), (b) At high bias, the trend is 
reversed and the electric field decreases with nanotube diam- 
eter. 
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FIG. 9: Potential profile of a nanotube at equilibrium but in 
the presence of an electric field corresponding to a bias of IV. 
Fermi energy is —0.5V throughout the tube. In contrast to 
the non equilibrium case shown in Fig|7] at equilibrium the 
potentials with and without scattering agree with each other. 
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FIG. 10: Current voltage characteristics of a (12,0) nanotube 
of different lengths. The dashed line shows the ballistic limit. 
The resistivity of the tube increases with the length causing 
current saturation. 
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FIG. 11: (a) Current-voltage characteristics of nanotubes of 
various diameters, in the presence of electron-phonon scatter- 
ing, (b) The differential conductance versus bias correspond- 
ing to (a). Inelastic phonon emission causes the decrease in 
differential conductance at high bias. Contribution of non 
crossing subbands leads to the increased low and high bias 
conductance for large diameter tubes. 
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FIG. 12: Schematics for tunneling into non crossing subbands 
at the source. The dashed line shows a potential profile V(y) 
and the dotted line shows a boundary between a contact and 
a nanotube. The arrow 1 shows injection from a perfect con- 
tact with a bandstructure of a nanotube. The arrow 2 shows 
injection from end-contacts with constant density of states or 
a scattering-assisted tunneling from crossing into non crossing 
subbands. 
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FIG. 13: The potential profile in a single- and multiwalled 
NT. The barrier width in the outer (240,0) shell decreases 
due to the additional screening provided by the inner shells. 
Inset: a comparison of differential conductance in singlewalled 
(240,0) nanotube (dashed) and multiwalled nanotube of the 
same diameter. Both curves are bell-shaped, signifying a 
quenching of Zener tunneling. 
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FIG. 14: Potential profile for good (a = 1) and poor (a = 
1CP 3 and a = 10 3 ) source contacts. When the coupling is 
poor, voltage drops mostly at the source, which opens Zener 
tunneling for higher subbands. 



0.14 
0.12 
— 0.1 
0.08 

CD 

™ 0.06 
° 0.04 
0.02 




L=213 nm 








0.4 



0.8 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8. 



a=1 

L=213 nm 







0.4 

V D (V) 



0.8 



FIG. 15: Differential conductance G vs. voltage Vd for a) 
poor (a = 10 -3 ) and b) good (a = 1) source contacts. Pres- 
ence of Zener tunneling at a = 10~ qualitatively changes 
G-V characteristics. 
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FIG. 16: Current as a function of a source contact quality 
factor a. When the contact is poor (far from 1), the effect 
of higher subbands becomes important which signifies the in- 
creasing trend of differential conductance versus bias. 



